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Abstract 

We study the dynamical phase diagram of a dilute Bose-Einstein conden- 
sate (BEC) trapped in a periodic potential. The dynamics is governed by a 
discrete non-linear Schrodinger equation: intrinsically localized excitations, 
including discrete solitons and breathers, can be created even if the BEC's 
interatomic potential is repulsive. Furthermore, we analyze the Anderson- 
Kasevich experiment [Science 282, 1686 (1998)], pointing out that mean field 

effects lead to a coherent destruction of the interwell Bloch oscillations. 
PACS: 63.20.Pw, 05.45.-a 
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Localization phenomena are ubiquitous in physics and in biology [0]. Intrinsically lo- 
calized excitations, as solitons (shape preserving) and breathers (characterized by internal 
oscillations), are important channels for energy transport in non-linear media, such as op- 
tical fibres and waveguides [Q, polaronic materials and biological molecules Intense 
theoretical research is now focusing on the existence of solitons and breathers in a lattice 
(often named discrete solitons and discrete breathers (DSB) (H), such as those appearing in 
quantum systems governed by a discrete non-linear Schrodinger equation (DNLSE) ||5|-0. 
Current approaches include the search for exact solutions in some limits P]; effective (point) 
particle and variational approaches |]9|-[TI|; perturbation around the linearized case and, of 
course, numerical solutions [0,0. Although intensely studied, DSB have been experimen- 
tally observed only quite recently in superconducting ladders of Josephson junctions 



m 



antiferromagnet systems [T^, in optical waveguides and in low dimensional materials 
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The discrete solitons/breathers are characterized by a dynamical, self-maintained energy 
localization, due to both the discreteness and the non-linearity of the underlying equations 
of motion. The discreteness provides a band structure of the excitation spectrum, while the 
non-linearity allows for the tuning of the DSB energy outside the band. The finite energy 
gap guarantees the (meta-) stability of the DSB's. These, obviously, have a different nature 
to the "Anderson localizations", created by impurities or imperfections of the lattice ||16|| ; 



the incorporation of disorder into non-linear excitations, and the tracing out of its dynamical 
effects are also important theoretical problems. 

Bright solitons can occur in spatially homogeneous, dilute Bose-Einstein condensates 
(BEC) with an attractive interatomic interaction (s-wave scattering length a < 0) p!7| , p!8 



Dark solitons, propagating density dips, have been predicted and experimentally observed 
in BEC's with a repulsive interaction (a > 0) |19[. The dynamics of a BEC trapped in a 



spatially periodic potential pO|-p4, on the other hand, can be mapped, in the tight binding 
approximation, to a DNLSE. Bright discrete solitons and breathers can be created in DNLSE 
even with a repulsive BEC interatomic interaction. 



In Ref. [^] a one-dimensional, vertical optical array was created by two counterprop- 
agating laser beams. A weakly interacting Bose-Einstein condensate was trapped in ~ 30 
wells, situated at the antinodes of the standing optical wave. Each well contained approx- 
imately 1000 condensate atoms, with the peak densities matching a gaussian profile. The 
lowest Bloch band dynamics of this system maps on a DNLSE, obtained by discretizing 
the Gross-Pitaevskii equation governing the condensate dynamics in the periodic potential. 
Since the array is oriented vertically, the atoms undergo coherent Bloch oscillations, driven 
by the interwell gravitational potential. At the edge of the Brillouin zone, a fraction of atoms 
can Zener tunnel in the higher energy band which, in this specific case, is in the continuum. 
A coherent leakage from the trap was observed at each Bloch period. However, an increase 
of the on-site particle densities washed out the coherent signal. This has been understood 



as a drift of the relative phases between wells due to mean field effects p3 . 

In this Letter we study the DNLSE using a variational approach: we show that a (com- 
plex) gaussian variational ansatz (describing a soliton / breather profile, combined with a 
Lagrangian optimization) yields a coupled dynamics of the profile parameters. A variety of 
discrete solitons / breathers are investigated analytically and compared with numerical cal- 
culations. This comparison is surprisingly successful in describing even details of the quite 
complex dynamical and collisional behaviour. Stability phase diagrams for such states are 
obtained by inspection of the profile dynamics equations. 

The full Bose condensate dynamics satisfies the Gross-Pitaevskii equation (GPE) p^ : 

= -|^v^*+[K.*+^oi$n$ (1) 

where Vext is the external potential and go = ^""^ with m the atomic mass. For a tilted 



trap as in |23| the external potential is given by the sum of the gravitational potential and 



the laser field 

Vext{r) = mgz + Ul{x, y) sin^[2'Kz / \] (2) 

where A is the wavelength of the lasers (the spacing in the lattice is A/2) and [/^(x, y) is 
determined by the transverse intensity profile of the (nearly gaussian) laser beams. In ||2^ 



A = 850 nm and the 1/e^ radius of the transverse profile is ~ 80/im, an order of magnitude 
larger than the transverse radius of the condensate. The well depths scale linearly with the 
intensity of the beam. At the center of the beam the trap depths are IAEr where Er = 
is the recoil energy {k = ^). 

In the tight-binding approximation the condensate order parameter can be written as: 

$(f,t) = ^Y.Mt) <P{r-rn), (3) 

n 

where is the total number of atoms and 0(r — r„) is the condensate wave function 



localized in the trap n with / dr (pn 4>n+i — 0, and J df (p"^ = 1. ipn = \J Pnif) e*^"*^*^ is the 
n-th amplitude {pn = Nn/Nx, where Nn and 6n are the number of particles and phases in 
the trap n). Replacing the ansatz in (|^) we find that the GPE reduces to a DNLSE: 

dl/jn 1 

i-Q^ = --{i^n-l + V'n+l) + {en + A\tpn H^n (4) 

where = ^ J df[^{V<Pn? + K.t^^], ^ = ^1 df<Pt and ~ - / df[^V<Pn ■ V</)„+i + 
0nKxt0n+i]; the time has been rescaled as t — ^> -^t. Eq. @ is the equation of motion 
tjjn = Q^^, ^ , where 7i is the Hamiltonian function 

with i?/)* , ipn canonically conjugate variables. Both the Hamiltonian Ti and the norm Yin I 
ipn P= 1 are conserved. 

To analyze the Anderson-Kasevich experiment, we study the dynamical evolution of 
a gaussian profile wave packet, which we parametrize as ipyif) = \fk ■ exp | — ^^Z^ — h 
ipin — ^) + — .^)^| where and 7(t) are, respectively, the center and the width 
(in the lattice units) of the density pn =| p(^) and 5(t) their associated momenta 

and A;(7,^) a normalization factor. The wave packet dynamical evolution can be obtained 
by a variational principle from the Lagrangian L = J2n ^'^nV'j^ — 'H, with the equations of 
motion for the variational parameters qi{t) = ^,7,^,5 given by = After some 

algebra we obtain C = pC, — ^ — + cosp ■ e'"^ — 1^(7, ^) where V = ^ + ^^'i 



V=k J dnenC ^ . The variational equations of motion become: 



^ = sinp ■ e ^ 



2A 4 dV 

(7) 



7 = 7^ cosp • e 

with the pairs ^, p and 5 which are canonically conjugate dynamical variables with respect 
to the effective Hamiltonian 

i7 = ^A_-cosp- e-^ + \/(7,0- (8) 

The wave packet group velocity is given hy Vg = ^ = ^ = tanp/m* with an inverse 
effective mass (m*)~^ = = cospe"''. The quasi-momentum dependence of the effective 
mass allows a rich variety of dynamical regimes. Solitonic solutions with a positive non- 
linear parameter A > 0, for instance, are allowed by a negative effective mass. A regime 
with a diverging effective mass m* oo leads to a self-trapping of the wave packet. 

In the following we will study the dynamical regimes of Eq. in two particular cases. 
We first consider a tilted (washboard) periodic potential describing the vertical optical trap 
created in the Anderson-Kasevich experiment. In particular, we will show that non-linear 
effects tend to destroy the Bloch oscillations, consistently with the experimental results. 



The external potential for the tilted trap is given by Eq. with e„ = oon where 



UJ 



mgX/2 
2K ■ 



We find V = uj^ and p = —oo . 



Then we consider the case of an horizontal array, in which gravitation only provides a 
constant energy shift, with p = 0. We will classify four different regimes which include 
discrete breathers and solitons. 

Tilted trap. It is well known that single atoms in a tilted washboard potential os- 
cillate among sites at the Bloch frequency. This regime is described by Eqs. |^) with 
A = (corresponding to a negligible mean field condensate interaction). This is, pre- 



cisely, the regime investigated in P3 in which a coherent output was observed. In- 



deed a variational estimate gives A ~ 0.5 and u; ~ 2 (the scaled time is in units of 



h/2K = 0.35ms). In this limit the equations of motion can be solved exactly giving 
p{t) = —ut + pQ, C,{t) = —A[cospo — cos{ujt — pq)], 7^(t) = AA'^logAuj{cos[2{ujt — pQ)] — 
cos2po — 4 sinpo sin (cijt — po) ^ 4sin^po} + 7o ^^id b(t) = — "^"^ lsmiojt-po)+smpo] ^]-^gj.g 
A = —Hq/uj cospo, 5(0) = 6o = and Hq is the (conserved) initial energy. In the inset of 
Fig. 1 we show that the Bloch oscillations described by the variational ansatz (solid line) 
are in excellent agreement with the full numerical solution of the DNLSE (dashed line). The 
numerical average position is defined as I]„n | 'ipn 

The effect of non-linearity on the Bloch oscillations is dramatic. This has been studied 
experimentally in by increasing the density in each well, and observing a degradation 
in the interference pattern. With A 7^ we have: 

The Eq. (||) displays an effective damping term proportional to the velocity ^. We stress that 
the dynamics is fully Hamiltonian, and real dissipative processes are absent. For t —>■ 00, 'j 
tends to a constant value 7/jn and 6 ~ ^\ t, so the term AS has the correct positive sign 

V 'fin 

(for large t). The apparent damping is the consequence of a diverging effective mass of the 



— i — ^ 

wave packet m* ~ e , which stops the Bloch oscillations. This effect manifests, in the 

numerical analysis of Eq. (^, as a distortion of the on-site phases. In Fig.l (solid line) we 
show the variational Bloch dynamics with a non-linear parameter A = 10 and initial values 



^0 = 0, = and 5q = 0. The oscillation roughly decreases as ^(t) ~ — A(l— e ^"^z*" cosLut). 
The dashed line shows the full numerical solution of the DNLSE, in good agreement with 
the analytical result (solid line). The discrepancy at t > 10 is due to the breaking of the 
gaussian wave packet in the numerical simulation. 

Untilted trap. The momentum p{t) = pq is conserved. We note that the equations of 
motion are invariant with respect to the replacement A —A, Po —>■ Po + n and t —t. 
In Fig. 2 we report the dynamical phase diagram as a function of A versus cosp, for 
A > and with initial values ^0 = and 60 = 0. This phase diagram has been calculated 
analytically from Eqs. (§-|l) and checked numerically. In the region cosp > there are two 
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distinct regimes. When i^o > 0, 7(t) < 'jmax- this is the self-trapped regime in which the 
boson wave packet remains locahzed around few sites. The self-locahzation is a genuine 
non-hnear effect, characterized by a diverging effective mass. In particular, the self-trapped 
wave packet cannot translate along the array: this is a major difference with respect to 
soliton-like solutions. The limit values for t ^ oo are 7 „ ^„ , 5 — 00 and f 0. We 
note that a non-linear self-trapping occurs also in a two-site problem |^ . 



A diffusive regime occurs when — cospo < -f^o ^ 0. In this case 7(t 00) 00 
and ^ ~ —Ho/ tanpo = const. The transition between the two regimes occurs at Ac = 
2^/^lJQ cospoe'^^^'^o _ With A > Ac, the ratio between the initial value of the width 70 and 
the limit width jmaxit —>■ 00) is given by 

7n A — Ac , . 

[max 

In Fig. 3 we plot Jo/'jrnax vs. A/Ac while in the inset we report the variational and numerical 
values of the width of the density vs. A/Ac after a time r = 10 (scaled units). The 
discrepancy at large A/Ac is due to a slight deviation of the numerical density profile from 
a gaussian shape. We checked the stability of self-trapping, also considering different initial 
forms of the wave packet. 

In the region cosp < the self-trapping condition is given by Hq >\ cospo I- iii this 
case Ac = 2-^7779 | cospo I (1 ~ e~'^^^'^o'^. A soliton solution can be determined by imposing 
j = S = 0. We find Asoi = 2^/77^^^ e'^/^'^o . For A = A^oi the center of the wave packet 
moves with a constant velocity ^ and its width remains constant in time. For Ac < A < Asoi, 
^ ^ 00 while 7(t) oscillates, corresponding to a breather solution. The numerical solution 
of Eq. d) confirms the variational predictions. In Fig. 4(a) we report the numerical density 
profile for A in the breather region with po = Sn/A; in Fig. 4(b) the numerical density 
profile is drawn for the same A and Pq = n/A: the change of the sign of cosp gives a breather 
solution in the first case, and the spreading of the density in the second one. In Figs. 4(c) 
and 4(d) we show that for A = A^o; the soliton after hitting a wall rebounds and regains its 
shape. 



In conclusion, we have studied the dynamics of a dilute Bose condensate trapped in a 
periodic potential. We have analyzed the Anderson-Kasevich experiment, pointing out that 
mean field effects lead to a coherent destruction of the interwell Bloch oscillations. We have 
shown that intrinsically localized excitations, like discrete breathers and solitons, can exist 
also with a repulsive interatomic potential by studying analytically their dynamical phase 
diagram. 

Discussions with S.R. Shenoy and P. Sodano are acknowledged. This work has been 
supported by the Cofinanziamento MURST. 
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FIGURES 

FIG. 1. Coherent destruction of Bloch oscillations: numerical (dashed line) and variational 
(solid line) average position of the density in the tilted trap (cj = 2, A = 10 and initial po = 0, 
6o = 0, 70 = 10). Inset: Bloch oscillations for A = 0. 

FIG. 2. Dynamical phase diagram in the untilted trap: A vs. cosp, for A > and with initial 
^0 = and 5o = 0. 

FIG. 3. Universal curve for Jo/jmax vs. A/A^ as in Eq. ([To|). The inset shows the variational 
(solid line) and numerical (circles) width of the density vs. A/Ac after the time r = 10, scaled 
units (initial values: po = 0, 5q = and 70 = 7). 

FIG. 4. (a) Numerical density profiles calculated at different times (t = 0,50,100,... ,400). 
The value of A is inside the breather region of the dynamical phase diagram 
{po = 37r/4, 7o = 10, A = 0.24); (b) the same parameters as in (a), but for pQ = 7r/4; (c) 
width (dotted line) and average position (solid line) calculated numerically for A = Agoi and 
Po = Stt/A in a finite array of 73 sites (the soliton hits a wall); (d) numerical density profile for 
t = 0, 50, 100, 150, 200 for the soliton described in (c). 
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